Search for: All records

Holmberg 15A (H15A), the brightest cluster galaxy of A85, has an exceptionally low central surface brightness even among local massive elliptical galaxies with distinct stellar cores, making it exceedingly challenging to obtain high-quality spectroscopy to detect a supermassive black hole (SMBH) at its center. Aided by the superb sensitivity and efficiency of KCWI at the Keck II Telescope, we have obtained spatially resolved stellar kinematics over a ∼100″ x 100″ contiguous field of H15A for this purpose. The velocity field exhibits a low-amplitude (∼20 km/s) rotation along a kinematic axis that is prominently misaligned from the photometric major axis, a strong indicator that H15A is triaxially shaped with unequal lengths for the three principal axes. Using 2500 observed kinematic constraints, we perform extensive calculations of stellar orbits with the triaxial Schwarzschild code, TriOS, and search over ~40,000 galaxy models to simultaneously determine the mass and intrinsic 3D shape parameters of H15A. We determine a ratio of p = 0.89 for the middle-to-long principal axes and q = 0.65 for the short-to-long principal axes. Our best estimate of the SMBH mass, M_BH = (2.16_{−0.18}^{+0.23})x10^{10} M⊙, makes H15A — along with NGC 4889 — the galaxy hosting the most massive SMBHs known in the local Universe. Both SMBHs lie significantly above the mean M_BH–σ scaling relation. Repeating the orbit modeling with the axisymmetrized version of TriOS produces worse fits to the KCWI kinematics and increases M_BH to (2.55 ± 0.20)x10^{10} M⊙, which is still significantly below the M_BH = (4.0 ± 0.8)x10^{10} M⊙ reported in a prior axisymmetric study of H15A. 

Abstract Supermassive black hole binaries (SMBHBs) are thought to form in galaxy mergers, possessing the potential to produce electromagnetic (EM) radiation as well as gravitational waves (GWs) detectable with pulsar timing arrays (PTAs). Once GWs from individually resolved SMBHBs are detected, the identification of the host galaxy will be a major challenge due to the ambiguity in possible EM signatures and the poor localization capability of PTAs. To aid EM observations in choosing follow-up sources, we use NANOGrav’s galaxy catalog to quantify the number of plausible hosts in both realistic and idealistic scenarios. We outline a host identification pipeline that injects a single-source GW signal into a simulated PTA data set, recovers the signal using production-level techniques, quantifies the localization region and number of galaxies contained therein, and finally imposes cuts on the galaxies using parameter estimates from the GW search. In an ideal case, the 90% credible areas span 29–241 deg², containing about 14–341 galaxies. After cuts, the number of galaxies remaining ranges from 22 at worst to one true host at best. In a realistic case, these areas range from 287 to 530 deg²and enclose about 285–1238 galaxies. After cuts, the number of galaxies is 397 at worst and 27 at best. While the signal-to-noise ratio is the primary determinant of the localization area of a given source, we find that the area is also influenced by the proximity to nearby pulsars on the sky and the binary chirp mass. 

Abstract Evidence for the majority of the supermassive black holes in the local Universe has been obtained dynamically from stellar motions with the Schwarzschild orbit superposition method. However, there have been only a handful of studies using simulated data to examine the ability of this method to reliably recover known input black hole massesM_BHand other galaxy parameters. Here, we conduct a comprehensive assessment of the reliability of the triaxial Schwarzschild method atsimultaneouslydeterminingM_BH, stellar mass-to-light ratioM*/L, dark matter mass, and three intrinsic triaxial shape parameters of simulated galaxies. For each of 25 rounds of mock observations using simulated stellar kinematics and theTriOScode, we derive best-fitting parameters and confidence intervals after a full search in the 6D parameter space with our likelihood-based model inference scheme. The two key mass parameters,M_BHandM*/L, are recovered within the 68% confidence interval, and other parameters are recovered between the 68% and 95% confidence intervals. The spatially varying velocity anisotropy of the stellar orbits is also well recovered. We explore whether the goodness-of-fit measure used for galaxy model selection in our pipeline is biased by variable complexity across the 6D parameter space. In our tests, adding a penalty term to the likelihood measure either makes little difference, or worsens the recovery in some cases. 

ABSTRACT We present wide-field, deep K-band photometry of 98 luminous early-type galaxies (ETGs) from the MASSIVE survey based on observations taken with the WIRCam instrument on the Canada–France–Hawaii Telescope. Using these images, we extract accurate total K-band luminosities (LK) and half-light radii (Re) for this sample of galaxies. We use these new values to explore the size–luminosity and Faber–Jackson relations for massive ETGs. Within this volume-limited sample, we find clear evidence for curvature in both relations, indicating that the most luminous galaxies tend to have larger sizes and smaller velocity dispersions than expected from a simple power-law fit to less luminous galaxies. Our measured relations are qualitatively consistent with the most massive elliptical galaxies forming largely through dissipationless mergers. When the sample is separated into fast and slow rotators, we find the slow rotators to exhibit similar changes in slope with increasing LK, suggesting that low-mass and high-mass slow rotators have different formation histories. The curvatures in the Re–LK and σ–LK relations cancel, leading to a relation between dynamical mass and luminosity that is well described by a single power law: Reσ2 ∝ LKb with b ≈ 1.2. This is consistent with the tilt of the fundamental plane observed in lower mass elliptical galaxies. 

Abstract The three-dimensional intrinsic shape of a galaxy and the mass of the central supermassive black hole provide key insight into the galaxy’s growth history over cosmic time. Standard assumptions of a spherical or axisymmetric shape can be simplistic and can bias the black hole mass inferred from the motions of stars within a galaxy. Here, we present spatially resolved stellar kinematics of M87 over a two-dimensional 250″ × 300″ contiguous field covering a radial range of 50 pc–12 kpc from integral-field spectroscopic observations at the Keck II Telescope. From about 5 kpc and outward, we detect a prominent 25 km s⁻¹rotational pattern, in which the kinematic axis (connecting the maximal receding and approaching velocities) is 40° misaligned with the photometric major axis of M87. The rotational amplitude and misalignment angle both decrease in the inner 5 kpc. Such misaligned and twisted velocity fields are a hallmark of triaxiality, indicating that M87 is not an axisymmetrically shaped galaxy. Triaxial Schwarzschild orbit modeling with more than 4000 observational constraints enabled us to determine simultaneously the shape and mass parameters. The models incorporate a radially declining profile for the stellar mass-to-light ratio suggested by stellar population studies. We find that M87 is strongly triaxial, with ratios ofp= 0.845 for the middle-to-long principal axes andq= 0.722 for the short-to-long principal axes, and determine the black hole mass to be $({5.37}_{-0.25}^{+0.37}\pm 0.22)\times {10}^9{M}_{\odot }$, where the second error indicates the systematic uncertainty associated with the distance to M87. 

Abstract We present a detailed analysis of the behavior of the triaxial Schwarzschild orbit superposition method near the axisymmetric limit. Orbit superposition modeling is the primary method used to determine dynamical masses of supermassive black holes ( M BH ) in nearby galaxies; however, prior studies have reported conflicting results when comparing the outcome from axisymmetric orbit codes with that from a triaxial orbit code in the axisymmetric limit. We show that in order to achieve (oblate) axisymmetry in a triaxial code, care needs to be taken to axisymmetrize the short-axis tube orbits and to exclude both the long-axis tube and box orbits from the orbit library. Using up to 12 Gauss–Hermite moments of the line-of-sight velocity distributions as constraints, we demonstrate the effects of orbit types on the best-fit M BH in orbit modeling of the massive elliptical galaxy NGC 1453 reported in Liepold et al. In addition, we verify the efficacy of our updated code on a mock galaxy data set. We identify a subset of slowly precessing quasi-planar orbits for which the typical integration times can be insufficient to fully capture the equilibrium orbital behavior in both axisymmetric and triaxial systems with central black holes. Further investigation is needed for a more reliable treatment of these orbits. 

Abstract Pulsar timing array observations have found evidence for an isotropic gravitational-wave background with the Hellings–Downs angular correlations between pulsar pairs. This interpretation hinges on the measured shape of the angular correlations, which is predominantly quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlations into a sum of Legendre polynomials and use a Bayesian analysis to constrain their coefficients with the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). When including Legendre polynomials with multipolesℓ≥ 2, we only find a significant signal in the quadrupole with an amplitude consistent with general relativity and nonzero at the ∼95% confidence level and a Bayes factor of 200. When we include multipolesℓ≤ 1, the Bayes factor evidence for quadrupole correlations decreases by more than an order of magnitude due to evidence for a monopolar signal at approximately 4 nHz, which has also been noted in previous analyses of the NANOGrav 15 yr data. Further work needs to be done in order to better characterize the properties of this monopolar signal and its effect on the evidence for quadrupolar angular correlations. 

Abstract We present a stellar dynamical mass measurement of a newly detected supermassive black hole (SMBH) at the center of the fast-rotating, massive elliptical galaxy NGC 2693 as part of the MASSIVE survey. We combine high signal-to-noise ratio integral field spectroscopy (IFS) from the Gemini Multi-Object Spectrograph with wide-field data from the Mitchell Spectrograph at McDonald Observatory to extract and model stellar kinematics of NGC 2693 from the central ∼150 pc out to ∼2.5 effective radii. Observations from Hubble Space Telescope WFC3 are used to determine the stellar light distribution. We perform fully triaxial Schwarzschild orbit modeling using the latest TriOS code and a Bayesian search in 6D galaxy model parameter space to determine NGC 2693's SMBH mass (M_BH), stellar mass-to-light ratio, dark matter content, and intrinsic shape. We find $M_{\mathrm{BH}}=\left(1.7\pm 0.4\right)\times {10}^9{M}_{\odot }$and a triaxial intrinsic shape with axis ratiosp=b/a= 0.902 ± 0.009 and $q=c/a={0.721}_{-0.010}^{+0.011}$, triaxiality parameterT= 0.39 ± 0.04. In comparison, the best-fit orbit model in the axisymmetric limit and (cylindrical) Jeans anisotropic model of NGC 2693 prefer $M_{\mathrm{BH}}=\left(2.4\pm 0.6\right)\times {10}^9{M}_{\odot }$and $M_{\mathrm{BH}}=\left(2.9\pm 0.3\right)\times {10}^9{M}_{\odot }$, respectively. Neither model can account for the non-axisymmetric stellar velocity features present in the IFS data.